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We present results of first-principles calculations of the frequency-dependent polarizabilities of all 
alkali atoms for light in the wavelength range 300-f 600 nm, with particular attention to wavelengths 
of common infrared lasers. We parameterize our results so that they can be extended accurately 
to arbitrary wavelengths above 800 nm. This work is motivated by recent experiments involving 
simultaneous optical trapping of two different alkali species. Our data can be used to predict 
the oscillation frequencies of optically-trapped atoms, and particularly the ratios of frequencies of 
different species held in the same trap. We identify wavelengths at which two different alkali atoms 
have the same oscillation frequency. 

PACS numbers: 31.15.Md, 32.10.Dk, 32.70.Jz, 32.80.Rm 
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I. INTRODUCTION 

The trapping of two different species of ultracold gases 
has been an active topic of research for about ten years 
(see pj and references therein), and has led to the simul- 
taneous production of quantum degenerate Bose-Einstein 
and Fermi-Dirac gases 0, 0, an d recently to the ob- 
servation of heteronuclear Feshbach resonances involving 
Bose-Einstein and Fermi-Dirac species 0, |(j . 

Two recent experiments 0, |(| provide a glimpse of 
the possibilities of mapping out the phase diagram of 
a Bosc-Einstcin/Fermi-Dirac mixture, potentially a rich 
system that remains to be explored. Our present work 
provides accurate estimates of the trapping frequencies of 
disparate species confined in the same optical trap. We 
have been attentive to finding cases in which all species, 
Bose-Einstein and Fermi-Dirac atoms, and the molecules 
formed from them, have the same trapping frequencies. 
This case is of particular interest because it has been 
shown to be optimal for a recently proposed scheme of 
cooling mediated by Feshbach resonances In addi- 
tion, we note that while the DC polarizabilities of alkali 
atoms have been subject to careful experimental inves- 
tigation, there is relatively little accurate experimental 
data on frequency-dependent polarizabilities. We draw 
attention to the fact that current experiments offer the 
possibility of measurements of unprecedented accuracy 
regarding the ratio of frequency-dependent polarizabili- 
ties of different species. Such measurements could pro- 
vide more decisive comparison of ab initio atomic theory 
with experiments, than has been made to date. 

We begin with a summary of the calculation of 
frequency-dependent polarizabilities, the details of which 
have been presented in previous publications 0, @- 
The computational method we use here has delivered 
agreement with experiments on static polarizabilities at 
the 1% level. There are few direct measurements of 



frequency-dependent polarizabilities - the experiments 
are much more demanding, and have larger uncertain- 
ties associated with modeling the intensity profile of a 
focused laser beam. However, the alkali metal atoms are 
a special case, at least in the infrared spectral region. 
There, the polarizability is dominated by the contribu- 
tion from the longest- wavelength resonance transitions, 
and it is possible to perform consistency checks between 
static values and long- wavelength values of the frequency- 
dependent atomic polarizabilities. We present tabulated 
values of calculated frequency-dependent polarizabilities 
over a range of wavelengths from the ultraviolet through 
the infrared spectrum, with particular attention to those 
wavelengths (primarily in the infrared) which have been 
or might be employed in the present generation of op- 
tical traps. For the infrared region, we display simple 
formulae that allow accurate computation of frequency- 
dependent polarizabilities for wavelengths not explicitly 
reported upon here. 

Finally, we enumerate selected wavelengths, at which 
two species of optically-trapped alkali atoms would have 
the same frequencies of oscillation in a common optical 
trap. Such wavelengths offer opportunities for mapping 
out phase diagrams of Bose-Einstein and Fermi-Dirac 
mixtures in a space of reduced dimensionality. Moreover, 
to the extent that molecular resonances consist of weakly 
bound states at large internuclear separation, the oscil- 
lation frequencies of trapped molecular species should be 
very close to the frequencies of their atomic components. 
This may present an ideal system for investigation. 



II. CALCULATION OF POLARIZABILITIES 

The frequency-dependent polarizability of an atom in 
its ground state may be separated into polarizability of 
the ionic core a C ore and the valence contribution a™ s . 
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TABLE I: The values of the frequency-dependent polarizabilities of alkali-metal atoms in ajj. The wavelengths (in air) are 
given in nm, the corresponding frequencies are given in a.u. 







Li 


Na 


K 


Rb 


Cs 


Ft 


1554 


0.02931 


201.0(7) 


189.4(2) 


381.2(8) 


424.0(7) 


570.9(1.2) 


407.9(2.0) 


1550 


0.02939 


201.3(7) 


189.6(2) 


381.8(8) 


424.7(7) 


572.2(1.2) 


408.5(2.0) 


1540 


0.02958 


201.9(7) 


190.0(2) 


383.4(8) 


426.6(7) 


575.6(1.2) 


410.1(2.0) 


1340 


0.03399 


218.1(8) 


201.0(2) 


428.2(9) 


479.5(7) 


673.7(1.3) 


455.0(2.1) 


1332 


0.03420 


219.0(8) 


201.6(2) 


430.7(9) 


482.5(7) 


679.5(1.3) 


457.5(2.1) 


1313 


0.03469 


221.1(8) 


203.0(2) 


436.9(9) 


489.9(7) 


694.1(1.3) 


463.7(2.1) 


1240 


0.03673 


231.0(8) 


209.3(2) 


465.9(9) 


524.6(7) 


764.6(1.4) 


492.8(2.2) 


1152 


0.03954 


247.0(8) 


219.4(2) 


516.1(1.0) 


585.7(8) 


899.9(1.6) 


543.6(2.3) 


1090 


0.04179 


262.7(8) 


228.8(2) 


568.7(1.1) 


650.8(8) 


1062.3(1.7) 


597.2(2.4) 


1064 


0.04281 


270.8(8) 


233.6(2) 


597.5(1.2) 


686.9(9) 


1162.1(1.9) 


627(3) 


1060 


0.04297 


272.1(8) 


234.3(2) 


602.4(1.2) 


693.1(9) 


1180.0(1.9) 


632(3) 


1053 


0.04326 


274.6(8) 


235.7(3) 


611.3(1.2) 


704.4(9) 


1213.2(1.9) 


641(3) 


1047 


0.04351 


276.7(8) 


237.0(3) 


619.3(1.2) 


714.5(9) 


1243.8(2.0) 


650(3) 


1030 


0.04422 


283.2(8) 


240.7(3) 


643.9(1.3) 


746.0(9) 


1343.6(2.1) 


675(3) 


985 


0.04624 


304.0(8) 


252.0(3) 


728.3(1.4) 


855.8(1.0) 


1759(3) 


765(3) 


980 


0.04648 


306.7(8) 


253.5(3) 


739.9(1.4) 


871.2(1.0) 


1828(3) 


778(3) 


y4o 


n ri/i qi k 


oz / .o^y ) 


Oft/1 /1/'Q\ 

zo4.4(o ) 


oof; a (~\ (\\ 

000.4(1.0 ) 


1UU1. 0(1.1 J 


Z0ol(4 ) 


OO0(0 ) 


930 


0.04898 


339.6(9) 


270.3(3) 


895.9(1.7) 


1084.4(1.2) 


3325(5) 


957(3) 


799.3 


0.05699 


549.7(1.1) 


353.8(4) 


3666(7) 


13483(14) 


-2333(4) 


-1474(11) 


700 


0.06507 


1983(3) 


553.1(0.6) 


-1394(3) 


-1192.7(1.4) 


-707.9(1.4) 


-3787(18) 


600 


0.07592 


-645.9(1.3) 


4506(5) 


-438.4(0.9) 


-421.8(0.7) 


-333.9(1.0) 


-523(3) 


500 


0.09110 


-200.1(0.8) 


-411.1(0.4) 


-201.3(0.5) 


-196.8(0.6) 


-166.3(1.0) 


-204(2) 


400 


0.11388 


-86.4(0.8) 


-134.0(0.2) 


-131.4(2.4) 


-109.7(1.1) 


-88.2(1.3) 


-111(4) 


300 


0.15183 


-38.4(1.1) 


-57.5(0.3) 


-47.2(0.4) 


-43.7(1.1) 


-35(2) 


-43(2) 



a.u., in which e/y/AiFeo, m c , and the reduced Planck con- 
stant h have the numerical value 1, in Eq.(J2J). Polar- 
izability in a.u. has the dimensions of volume, and its 
numerical values presented here are thus measured in 
units of afj, where eto ~ 0.052918 nm is the Bohr ra- 
dius. The atomic units for a can be be converted to SI 
units via a/h [Hz/(V/m) 2 ]=2.48832x 10" 8 a [a.u.], where 
the conversion coefficient is ^Tre^a^/h and Planck con- 
stant h is factored out. The atomic unit of frequency u 
is E h /h « 4.1341 x 10 16 Hz, where E h is the Hartree 
energy. 

The sum over n! in Eq.Q converges rapidly, and un- 
less the frequency is resonant with the transitions other 
than the primary ones (ns — npi/2 and ns — np 3 / 2 )i the 
first term is dominant. Therefore, only a few terms need 
to be calculated accurately We use the experimental 
numbers compiled in Ref. 1 1 1 1| together with their uncer- 
tainties for the first term, for example for n! = 3 term for 
Na. High-precision theoretical all-order values are used 
for the next three terms (for example n' = 4, 5, 6 terms 
for Na). We refer the reader to Refs. [HE El for the 
detailed description of the single-double (SD) all-order 
method and its extensions. The all-order values of the 
matrix elements are evaluated for their accuracy, and the 
extensions of the SD all-order method (SDpT, which par- 
tially includes triple excitations or SD SC , which includes 
semi-empirical scaling of dominant terms) are used for 
certain transitions where those values are expected to 
be of better accuracy. The uncertainty of the resulting 
matrix elements is evaluated based on the relative impor- 



The calculation of the core polarizability is carried out in 
the random-phase approximation (RPA); we verified that 
our static values agree with RPA values from Ref. [TfT |. 
The separation of the polarizability into core and valence 
parts also requires the addition of the compensation term 
a vc which accounts for the contribution from the excita- 
tion to the occupied valence shell that is forbidden by 
the Pauli exclusion principle. The calculation of the core 
polarizability does not exclude the excitation from the 
core to the valence shell and half of this contribution has 
to be subtracted. Thus, the polarizability contributions 
may be separated as 

a ns = a coie + a vc + a™. (1) 

We calculate the valence part of the ac polarizability of 
the ground state ns for the alkali-metal atoms by com- 
puting the sum over states 



Ins) 2 



' SE n/pi/2 (n'p 1/2 \\D\ 
3 ^ \ SE 2 lri - lj 2 

n' \ n Pi/2 

SE n , p3/2 (n'p 3/2 \\D\\ns) 2 \ 

" P3/2 / 

where (np\\D\\ns) is the reduced elcctric-dipole matrix 
element, 5E n , pl/2 = E n > Pl/2 - E ns , and 5E n > P3/2 = 
E n ' P3/2 — E ns . In this formula, w is assumed to be 
at least several linewidths off resonance with the corre- 
sponding transition. We use the system of atomic units, 
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TABLE II: The parameters A na (in a.u.) for the calculation of ground state frequency-dependent polarizabilities of alkali- 
metal atoms in a 3 ) for the wavelengths above 800 nm. The corresponding uncertainties in the values of A" s are given in the 
last column. The corresponding reduced matrix elements and energy differences are given in columns labeled (np||D||ns) and 



Atom 


n 


(npi/2 \\D\\ ns) 


(nps/2 \\D\\ns) 


5E, Wl/2 






SA na 


Li 


2 


3.317(4) 


4.689(5) 


0.067906 


0.067907 


2.04 


0.69 


Na 


3 


3.5246(23) 


4.9838(34) 


0.077258 


0.077336 


1.86 


0.12 


K 


4 


4.102(5) 


5.800(8) 


0.059165 


0.059428 


6.26 


0.33 


Rb 


5 


4.231(3) 


5.977(4) 


0.057314 


0.058396 


10.54 


0.60 


Cs 


6 


4.4890(65) 


6.3238(73) 


0.050932 


0.053456 


17.35 


1.00 


Fr 


7 


4.277(8) 


5.898(15) 


0.055758 


0.063442 


24.8 


1.8 



tance of certain classes of the all-order terms, extensive 
comparison of various atom ic p roperties of alkali-metal 
atoms with experiment 0, ^[ 0], and the spread of 
SD, SDpT, and SD SC values. Some of these values were 
published previously in Refs. [TTL Il4j. The experimental 
values from are used for 6s — 7pi/ 2 and 6s — 7p 3 / 2 
transitions. In summary, the set of values used for the 
calculation of the first four terms in the sum over states 
of Eq.(J2J) consists of the best known values for these tran- 
sitions. The experimental energies from Refs. [l5llT^| are 
used for the corresponding energy levels. We refer to the 
total contribution from the first four terms as the main 
term; the remaining contributions are referred to below 
as the tail contribution; i.e. 



,(n' = 



+ 3) +a tal i(n' > ri + 3). (3) 



We evaluated the tail contributions with a complete 
set of basis Dirac-Hartree-Fock (DHF ) w ave functions 
generated using the B-spline method [13] • We use 40 
splines of order k — 7 for each angular momentum. The 
basis set orbitals are defined on a non-linear grid and 
are constrained to a large spherical cavity of a radius 
R = 75—110 a.u. The cavity radius is chosen to ac- 
commodate first four ns, np 1 / 2 , and np 3 / 2 valence or- 
bitals. For consistency, we use the same basis sets for 
the tail calculations as we used in the all-order calcula- 
tion of the main term matrix elements. The tail con- 
tribution is calculated in the DHF approximation. Wc 
list our results for the frequency-dependent polarizabil- 
ities from the ultraviolet through the infrared spectrum 
in Table |U Particular attention is given to the values 
for the infrared wavelengths which have been or might 
be employed in the present generation of optical traps. 
Our values of the ac polarizabilities are accurate to better 
than 0.5% with the exception of the values correspond- 
ing to wavelengths in the ultraviolet spectrum. For these 
wavelengths, non-primary resonances become dominant 
and the uncertainty of the polarizability values is domi- 
nated by the uncertainty in the corresponding resonance 
matrix elements. For example, the first resonances in Rb 
(5s — 5p transitions) occur at 795 nm and 780 nm, while 
next resonances (5s — 6p) transitions are near 420 nm. As 
the accuracy of the 5s — 6p matrix elements is lower than 
that of the 5s — 5p matrix elements the resulting accuracy 
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FIG. 1: The ratios a(A)/m of frequency-dependent polariz- 
abilities to the atomic weig hts for 6 Li and 23 Na. 

of the ac polarizabilities at 400 nm is also lower. 

III. RESULTS 

We have investigated the dependence on the frequency 
u! of all of the terms contributing to a ns , where ns is the 
ground state. We find that the a cor e, ct vc , atail, and a" s 
terms with n' = n+l,n + 2,n + 3 depend weakly on the 
value of u> for wavelengths above 800 nm. In fact, it is 
possible to accurately parameterize the values of the ac 
polarizabilities for those wavelengths by the formula 



l 8E npi/2 {np 1/2 \\D\\nsf 

l 5E np 3/2 (np 3/2 \\D\\nsy 
3 



(4) 



where A ns is independent on the value of ui. All val- 
ues in Eq.QJ are in atomic units. The values of the 
constants A ns and their uncertainties are given in Ta- 
ble HTl The values of A ns are calculated as the aver- 
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TABLE III: Selected wavelengths A a ; r in nm, at which two 
species of optically-trapped alkali atoms ( 6 Li, 7 Li, 23 Na, 40 K, 
41 K, 87 Rb, and 133 Cs) would have the same frequencies of os- 
cillation in a common optical trap. The corresponding ground 
state frequency-dependent polarizabilities of the alkali-metal 



atoms (in a.u.) 


are also 


given. 








1 




2 




557.56(13) 


b Li 


-360(1) 


2a Na 


-1375(1) 


768.648(3) 


6 Li 


681.3(1.3) 


40 R 


4525(112) 


799.46(12) 


6 Li 


549(1) 


40 K 


3649(7) 


768.645(3) 


6 Li 


681.3(1.3) 


41 K 


4638(112) 


798.44(12) 


6 Li 


553(1) 


41 K 


3763(7) 


786.07(1) 


6 Li 


597.9(1.2) 


87 Rb 


8639(20) 


804.57(4) 


6 Li 


533(1) 


87 Rb 


7704(8) 


862.097(28) 


6 Li 


412.5(1.0) 


133 Cs 


9113(22) 


904.76(4) 


6 Li 


361.7(9) 


133 Cs 


7991(12) 


549.72(18) 


7 Li 


-329(1) 


23 Na 


-1078(1) 


768 665 f 3") 


7 Li 


681(1) 


40 K 


3880(1 12") 


806.9(2) 


7 Li 


526(1) 


40 K 


2998(6) 


768 669 CX\ 


7 Li 


681(1) 


41 K 


3978n 1 9\ 


805.55(15) 


7 Li 


530(1) 


41 K 


3096(6) 


i oO.OO^ 1 ) 


7 T i 

J_j1 


oyo.y ^i.z j 


87 Rh 


( oooyiij ) 


807 16(4^ 


7 Li 


526(1) 


87 Rb 






.LI 


41U.0^1.U J 


133 p c 
OS 


7777 '(Of}} 
lilt \Z\J ) 


Q07 9^6^ 


7 Li 


359 3('Q , I 

OiJc/ .O Iff 


133 Cs 


UOUU 1 1U 1 


70Q OA(i(7\ 


23 NT a 
In a 


O04.y^4 ) 


87 Rh 




y4u.oo4^o ) 


23 Nn 

in a 


Z04. 1 ^OJ 


87 Rh 
xxu 


QQQM ^ 

yyy^i ) 


875.53(3) 


IN a 


295.5(3) 


l33r->„ 


1709(11) 


1091 A\ 
1UZ 1.014 ) 


23 Na 

in a 


9/19 Pifl} 


133p 




768.848(3) 


40 t,' 
K 


-3890(110) 


Kb 


-8470(10) 


78/1 qoM^ 

i o4.oZ^ 1 J 


40 K 


f^snftfi ^ 

DoUD^lO J 


87 Rh 

TV IJ 


1 A$IV\V\(97\ 
14oUU^Z / ) 


807 *31 ( (C\ 


l\ 


zyuol^O ) 


87nu 
IX D 


0404^0 ) 


868.90(4) 


40 K 


1298(3) 


133 Cs 


4317(14) 


938.6(2) 


40 K 


862(2) 


133 Cs 


2868(4) 


768.851(3) 


41 K 


-4013(116) 


87 Rb 


-8472(10) 


784.42(1) 


41 K 


6767(12) 


87 Rb 


14358(26) 


808.07(7) 


41 K 


2916(6) 


87 Rb 


6187(6) 


869.14(4) 


41 K 


1296(2) 


133 Cs 


4203(14) 


940.9(2) 


41 K 


854(2) 


133 Cs 


2771(4) 


790.303(6) 


87 Rb 


-1299(17) 


133 Cs 


-1986(3) 


873.39(4) 


87 Rb 


1624(2) 


133 Cs 


2483(12) 


1158(3) 


87 Rb 


580.5(8) 


133 Cs 


888(2) 



age values of the totals of the a core , a vc , a>t a u, and 
a™ s terms with n' = n + l,n + 2,n + 3 calculated in 
the range from 800 nm to 1600 nm. The values of the 
(np\\D\\ns) reduced electric-dipole matrix elements and 
the corresponding energy differences 5E are also given. 
The electric-dipole matrix elements for Li are experimen- 
tal values from Ref. an other matrix elements are 
experimental values compiled in [ill. The energy differ- 
ences are taken from Refs. 0, llfij- The values of the ac 
polarizabilities obtained using Eq.Q differ from our re- 
sults in Table [I] in the relevant wavelength range by less 



800 
600 • 
400 • 
200 


■200 
■400 
■600 



-800 



768.838nm 



740 




87 Rb 
40 K 



820 



FIG. 2: The ratios a(A)/m of frequency-dependent polariz- 
abilities to the atomic weig hts for 40 K and 87 Rb. 



than 0.1% which is significantly below the accuracy of 
the values themselves. We note that the formula |0J is 
not applicable for the lower wavelengths owing to the im- 
portance of the other n'p — ns resonances. For example, 
the 420 nm wavelength corresponds to 6p — 5s transition 
in Rb; therefore, the ground state polarizability at this 
frequency will be dominated by the n' = 6 term. 



IV. THE FREQUENCY-MATCHING CRITERIA 
FOR MIXED SPECIES 



We have also conducted a search of the wavelength val- 
ues at which two species of optically-trapped alkali atoms 
( 6 Li, 7 Li, 23 Na, 40 K, 41 K, 87 Rb, 133 Cs, and 210 Fr) would 
have the same frequencies of oscillation in a common op- 
tical trap, i.e. A such that 



'ai(A) 11J2 
mi a 2 (A) 



1. 



(5) 



where m; is the atomic weight |T^. The wavelengths of 
those selected matches and their uncertainties are listed 
in Table IlTTI We have omitted all matches at wavelengths 
below 500 nm. While we had found numerous matches 
in the range from 300 nm to 500 nm they are difficult 
to place accurately since the accuracy of our calcula- 
tion is lower in this region. The width of the matches 
in the 300 — 500 nm region is generally very narrow and 
the corresponding values of the polarizabilities are small, 
making them of limited experimental use. We define the 
width of the match as the wavelength range where the 
parameter s given by Eq. © is between 0.9 and 1.1. The 
width of the match is very narrow if the corresponding 
wavelength is in the immediate vicinity of the resonance. 
We have also omitted most matches at higher wave- 
lengths from Table II I II that are close to the resonance. 
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Note that we listed the values of the wavelengths in the 
air. While the difference between air and vacuum wave- 
lengths is small in the considered region, 0.15 — 0.4 nm, it 
is sufficient to cause changes in the frequency-dependent 
polarizabilities (at some wavelengths) which exceed the 
uncertainties of the calculations. We used the following 
formula to compute the wavelength in air from the vac- 
uum wavelength [20j 



A a 



A„/(l + 10" 8 (8060.51 



2480990 



132.274 



10 6 



17455.7 



39.32957 -a 2 x 10 6 



))• 



(6) 



The wavenumber a is here to be expressed in reciprocal 
nm. We also list the values of the ac polarizabilities and 
their uncertainties in Table ITTT1 

The uncertainty in the wavelength of the match results 
from the uncertainties in the values of the frequency- 
dependent polarizabilities ai(A) and «2(A) given in Ta- 
ble lHIl The uncertainties in the values of a(A) are found 
to be independent of A within the range of the match to 
a very good precision but are generally different for a\ 
and a-i- We calculate the resulting uncertainties <5Ai and 
8\i due to uncertainties in the values of ot\ and a.2 sep- 
arately and add the results in quadrature to obtain the 
final uncertainty 8X. To evaluate <5Ai, we calculate the 
a-y{ui)/mi = [c*i(A) ± 5a\]/mi and then determine the 
wavelengths A ± at which af(X)/mi = a2{X)/m2- The 
maximum of the differences between the initial match 
wavelength and A is taken to be the uncertainty 5Xi. 
The evaluation of 6X2 is done in the same way. 

We illustrate several of the matches given in Tabic [ffi] 
in Fig. ^ an d Fig- H The ratios a(X)/m of frequency- 
dependent polarizabilities to the atomic weights for Li 
and Na in the wavelength range 500-800 nm are plotted 
in Fig. ^ The corresponding matches are shown by the 
arrows and the corresponding wavelengths are given on 



the plot. The values ^ for 6 Li and 23 Na match at three 
wavelengths in the wavelength range plotted on Fig. ^ 
Two of the matches, at 589.397 nm and 670.79 nm, are 
very close to the resonances and, therefore, are not listed 
in Table IIIII as discussed before. The match areas for 
those cases are extremely narrow and the polarizability 
of the corresponding alkali is very sensitive to very small 
changes in the wavelength. Fig. |3 shows similar plot to 
illustrate the wavelengths at which 40 K and 87 Rb have 
the same frequencies of oscillation in a common optical 
trap in the range 740-820 nm. 



V. CONCLUSIONS 

In summary, we conducted a systematic study of the 
ground state frequency-dependent polarizabilities of the 
alkali atoms from ultraviolet through infrared spectral 
regions. The values of the ac polarizabilities and their 
uncertainties are calculated for a number of wavelengths, 
including the wavelengths of common infrared lasers. 
A combination of high-precision measurements of the 
ratios of frequency-dependent polarizabilities of different 
species could provide excellent tests of the current 
experimental and theoretical values of the elcctric-dipole 
matrix elements in alkali-metal atoms. We provide 
formulas and the necessary parameters for the accurate 
calculation of the ac polarizabilities for all alkali-metal 
atoms at wavelengths above 800 nm. Finally, we list 
selected wavelengths at which two species of alkali atoms 
would have the same oscillation frequencies in a common 
optical trap. 
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